改进的同步迭代算法在光声血管成像中的应用

光声成像结合了光学成像和声学成像的优点，是一种高分辨率，高对比度的无损伤医学成像技术.一种改进的同步迭代算法应用于光声图像重建.仿真和模拟结果表明，与传统的代数迭代算法相比，在90°，135°，180°的有限场光声成像中，此算法对测量误差的校正和迭代次数的收敛上具有较大的优势，图像重建的速度和成像质量都有了明显的提高.实验中，一种圆形扫描结构的光声成像装置，用于180°的有限场扫描，利用改进的同步迭代算法，重建出了高对比度和高分辨率（60μm）的鸡胚胎光声血管图像.实验证明，这种算法的应用，大幅度减少了数据采集时间，为光声成像技术运用于实时监测血流灌注和肿瘤光动力治疗的血管损伤效应提供了潜在的应用价值. 关键词： 光声成像 有限角度 代数迭代算法 光声血管成像     

用三角函数法获得非线性Boussinesq方程的广义孤子解

找到一个合适的代换——三角函数法，将非线性Boussinesq微分方程转换为非线性代数方程组.用吴消元法求解该非线性代数方程组，从而获得一般形式Boussinesq微分方程的广义孤子解. 关键词： Boussinesq方程 吴消元法 非线性代数方程组 孤子解     

可积开边界条件下的q形变玻色子模型

利用代数Bethe Ansatz方法在可积开边界条件下推广了q形变玻色子模型，得到可积开边界条件下此模型的哈密顿量及其本征方程.该工作可为在更小尺度下研究具有相互作用的玻色子系统提供有效的理论基础. 关键词： 代数Bethe Ansatz q形变玻色子模型')" href="#">q形变玻色子模型 开边界 可积系统     

Implicit meanflow–multigrid algorithms for Reynolds stress model computation of 3‐D anisotropy‐driven and compressible flows

The present paper investigates the multigrid (MG) acceleration of compressible Reynolds‐averaged Navier–Stokes computations using Reynolds‐stress model 7‐equation turbulence closures, as well as lower‐level 2‐equation models. The basic single‐grid SG algorithm combines upwind‐biased discretization with a subiterative local‐dual‐time‐stepping time‐integration procedure. MG acceleration, using characteristic MG restriction and prolongation operators, is applied on meanflow variables only (MF–MG), turbulence variables being simply injected onto coarser grids. A previously developed non‐time‐consistent (for steady flows) full‐approximation‐multigrid (s–MG) is assessed for 3‐D anisotropy‐driven and/or separated flows, which are dominated by the convergence of turbulence variables. Even for these difficult test cases CPU‐speed‐ups r_CPUSUP∈[3, 5] are obtained. Alternative, potentially time‐consistent approaches (unsteady u–MG), where MG acceleration is applied at each subiteration, are also examined, using different subiterative strategies, MG cycles, and turbulence models. For 2‐D shock wave/turbulent boundary layer interaction, the fastest s–MG approach, with a V(2, 0) sawtooth cycle, systematically yields CPU‐speed‐ups of 5±½, quasi‐independent of the particular turbulence closure used. Copyright © 2008 John Wiley & Sons, Ltd.     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

A new general algorithmic method in reaction syntheses using linear algebra

We discuss here a new general linear algebraic method (both model and algorithm) for describing and generating (among others) minimal reactions and also minimal mechanisms in stoichiometry, or dimensionless groups in physics as well. (Further applications in process network syntheses will be discussed in .) With some minor modifications of the input this method can be extended for several related questions: for generating direct and overall reactions, direct (steady state) mechanisms, for finding the possible resulting (overall) reactions among all possible mechanisms, etc.Computational results in section 4 show the speed of our algorithm.We give also mathematical background and results in sections 3, 5 and 6. However, we do not restrict ourselves to mathematics only, we also talk on the language of chemistry, too.The theoretical results in sections 3.2, 3.3, 5 and the computational examples in section 4 are completely new, further theoretical results will appear in and in .     

A multigrid procedure for Cartesian ghost‐cell methods

This paper proposes a multigrid technique for Cartesian grid flow solvers. A recently developed ghost body‐cell method for inviscid flows is combined with a nested‐level local refinement procedure, which employs multigrid to accelerate convergence to steady state. Different from standard multigrid applications for body‐fitted grids, a fictitious residual needs to be defined in the ghost cells to perform a correct residual collection and thus to avoid possible stalling of the multigrid procedure. The efficiency of the proposed local refinement multigrid Cartesian method is demonstrated for the case of the inviscid subsonic flow past a circular body. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180‐degree turn by algebraic Reynolds stress model

Turbulent flow in a rectangular duct with a sharp 180‐degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure‐driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large‐scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180‐degree turn using the algebraic Reynolds stress model. A boundary‐fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser‐Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log‐law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180‐degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log‐law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

抛物问题的mortar有限元逼近的瀑布型多重网格法

用瀑布型多重网格法解决椭圆、抛物问题，已有不少研究工作[1-2]，本文对抛物问题的mortar有限元的全离散格式提出瀑布型多重网格法，证明了该方法是最优的，即具有最优精确度和复杂度．     

Local cyclic homology

Local (cyclic) homology of filtered associative algebras is defined and discussed. In the commutative case, it fits naturally in the setting of (formal) de Rham cohomology; we thus obtain a satisfactory complement to a result of Feigin and Tsygan. There is a spectral sequence from discrete (cyclic) homology of the associated graded algebra to local (cyclic) homology of the original filtered algebra; we exhibit manageable conditions on its convergence.